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Abstract 

Current experimental data on the neutrino parameters is in good agreement with 
tribimaximal mixing and may indicate the presence of an underlying family symmetry. 
For 76 flavor groups, we perform a systematic scan for models: The particle content is 
that of the Standard Model plus up to three flavon fields, and the effective Lagrangian 
contains all terms of mass dimension < 6. We find that 44 groups can accommodate 
models that are consistent with experiment at 3a, and 38 groups can have models 
that are tribimaximal. For one particular group, we look at correlations between 
the mixing angles and make a prediction for #13 that will be testable in the near 
future. We present the details of a model with 12 = 33.9°, 6» 23 = 40.9°, 9 13 = 5.1° to 
show that the recent tentative hints of a non-zero 613 can easily be accommodated. 
The smallest group for which we find tribimaximal mixing is Tj. We argue that T7 
and T13 are as suited to produce tribimaximal mixing as A^ and should therefore be 
considered on equal footing. In the appendices, we present some new mathematical 
methods and results that may prove useful for future model building efforts. 
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1. Introduction 



Neutrino physics is a fast developing field. The past decade has seen the discovery of 
neutrino masses [1,2] and ever improving measurements of the neutrino mixing matrix 
Upmtsss [3,4]. Our growing knowledge of the neutrino parameters [5,6] has almost raised 
more questions than it answered: Why are neutrinos so light? Why are two of the mixing 
angles large and one vanishingly small? Why is C/pmns so different from Uckm [7]? These 
are some of the questions that any model for the neutrino sector needs to address. 
Experimental data suggests that the mixing angles are in good agreement with tribi- 
maximal mixing (TBM) [8,9]. The very form of the Harrison-Perkins-Scott matrix £7hps 
is suggestive of an underlying family symmetry between the three generations of leptons. 
In the past years, much effort has been vested in finding a family symmetry that would 
naturally lead to tribimaximal mixing, and to that end, some twenty odd groups have 
been the subject of model building efforts (see refs. [10-12] and references therein). 

It has been argued that A4 is particularly relevant for producing tribimaximal mixing 
[13-16], and by the number of publications (see e.g. Tab. 2 in ref. [10]) it is certainly the 
most popular discrete symmetry used for model building. That is why we start out by 
following down the same path to construct all A4 x Z3 models with up to three flavon 
fields where the lepton doublet L transforms as a triplet. We find 22,932 inequivalent 
models, of which 4,481 (19.5%) give mixing angles that are consistent with experiment at 
3a, and 4,233 (18.5%) that are tribimaximal. Restricting # 12 and 823 to their respective 
3a intervals, we obtain an interesting prediction for # 13 whose value is currently not 
known with very high precision: The by far most likely value is #13 = 0°, and there are 
extremely few models for 0° ^ #13 < 12°. We also present a model where all three mixing 
angles #12 — 34°, 623 — 41° and #13 ~ 5° lie in their respective la intervals to show that 
it is possible to accommodate the recent tentative hints of a non-zero #13 [6, 17, 18]. 

Next we explore whether is really special or we are looking for tribimaximal models 
"under the lamppost". There are 1,048 groups with less than or equal to 100 elements, 
and 90 of them have a 3-dimensional irreducible representation. For 76 groups, we 
construct all models with up to three flavon fields where the lepton doublet L transforms 
in a 3-dimensional irreducible representation. For the remaining 14 groups, a systematic 
scan would simply take too long. We find 44 groups (58%) that can accommodate 
models which are consistent with experiment at the 3a level, and 38 groups (50%) that 
can produce tribimaximal mixing. The smallest group for which we find tribimaximal 
mixing is T7, and the group with the largest fraction of tribimaximal models is T13. 
Incidentally, for T13 (and the other metacyclic groups) the set of tribimaximal models 
and the set of 3a models are almost identical, and this may be pointing towards a 
profound connection between T 13 and tribimaximal mixing that is more pronounced as 
compared to A4. For a recent publication that uses T 13 for model building, see ref. [19]. 

For our analysis, the computer algebra program GAP [20] played a central role. We 
used GAP to obtain the character table, the dimension of the conjugacy classes and the 
explicit form of the representation matrices for the 76 groups that we considered in 
this publication. In contrast to e.g. solving renormalization group equations, the use 
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of computers for algebraic and group theoretic operations is not widespread (a notable 
exception is ref. [21]). We strongly advocate the use of the SmallGroups Library [22] 
which collects in one place and provides easy access to all finite groups of order at most 
2,000 (except 1,024). 

In the appendices we present some new developments and mathematical background 
information relevant for model building with discrete symmetries. 

In Appendix A, we list the 90 groups of order less than or equal to 100 that have a 
3-dimensional irreducible representation. For each group, we indicate whether it is a 
subset of U(3), U(2) or U(2) x U(l), and at the same time check whether it contains A 4 
as a subgroup. Due to its length, the full list of the 1,048 groups of order at most 100 
is presented in a separate file [23]. 

In Appendix B we give the full details on how we generated the 1,048 groups and 
compiled the tables in Appendix A and ref. [23]. We elaborate on some disagreement 
that we have with the existing literature. 

In Appendix C we show how to find the vacuum expectation values that break a given 
group to any one of its subgroups. Unfortunately, finding all possible symmetry breaking 
patterns does not allow us to classify the models, since different vevs inducing the same 
symmetry breaking chain may lead to different mixing angles. 

In Appendix D, we discuss an algorithm due to van den Broek and Cornwell [24] for cal- 
culating the Clebsch-Gordan coefficients for any finite group. This allows us to construct 
the group invariants, or more generally, contract the family indices in the Lagrangian 
without referring to heuristic constructions as is common practice in the current litera- 
ture. 

Finally, in Appendix E we outline some of the most important concepts and theorems 
from the theory of groups that pertain to the present publication. 

2. Experimental Constraints 

The leptonic mixing matrix ?7pmns is generally parametrized by three angles, $i 2 , 9 23 , 
#13, and one Dirac phase S [25]. If the neutrinos are Majorana particles, there are two 
extra phases (f)% and $2 that do not affect neutrino oscillation phenomena [26] and are 
likely to remain unconstrained in the near future. In this paper, we use the standard 
parametrization [25] of £7pmns except for the definition of the Majorana phases, where 
we follow ref. [27]: 

/ C12C13 S12C13 si 3 e l5 \ 

C^pmns = -S12C23 - ci 2 si3S23e i<5 ci 2 c 2 3 - si2Si3S23e 4<5 C13S23 • diag^* 1 , e^ 2 , 1) (1) 

\ S12S23 - C 12 Si 3 C 23 e tS -C12S23 - Sl2Sl3C23e i<5 C13C23 / 

For comparing our results from Section 5 to experiment, we used refs. [5,6,17]. In Tab. 1, 
we summarize the relevant information for the reader's convenience. 
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Parameter 


Mean Value 


lcr range 


3a range 


#12 


34.4° 


33.4° - 35.4° 


31.5° -37.6° 


#23 


42.8° 


39.9° -47.5° 


35.5° - 53.5° 


#13 


5.6° 


2.9° -8.6° 


0° - 12.5° 



Table 1: The leptonic mixing angles from the global fit to data from ref. [6] (first table, left column). 

The solar and atmospheric neutrino mixing angles, #12 and #23, are relatively well- 
determined. #13, on the other hand, effectively only has an upper bound. The ex- 
perimental data is consistent with #13 being zero, as e.g. in exact tribimaximal mixing. 
If #13 = 0°, the Dirac phase loses physical significance. Currently, there are possible 
hints for a non-zero #13 [6, 17, 18]. A new generation of neutrino experiments will probe 
sin 2 #13 down to about 1CT 2 [17]. 

3. A Paradigm: A4 x C3 Family Symmetry with Three 
Flavon Fields 

To illustrate our general approach, we will choose A4 x C3 as the family symmetry and 
reproduce the results of the now classic paper by Altarelli and Feruglio [16]. Here and in 
the following we will use the alternate notation C n for Z n . Note that we could have taken 
any of the 439,820 models that we will be constructing in Section 4, but considering a 
model that is already well-known has the advantage of a clearer presentation of our 
methodology by stressing the differences to other approaches. 

The following lines of GAP code give us information on the group A 4 x C3: 

1 group := SmallGroup(36 , 11) ; ; 

2 Display (StructureDescription(group) ) ; 

3 chartab := Irr (group) ; ; 

4 Display (chartab) ; 

5 SizesConjugacyClasses (CharacterTable (group) ) ; 

6 LoadPackageO'repsn") ; ; 

7 for i in [1 . .Size(chartab)] do 

8 rep := IrreducibleAf f ordingRepresentat ion (chartab [i] ) ; 

9 for el in Elements (group) do 

10 Display (el'rep) ; 

11 od; 

12 od; 

These lines can be entered directly at the GAP prompt or saved in a file and executed 
automatically as explained later. Line 1 defines the group in terms of its GAPID (see 
Section B.l). Lines 4 and 5 display the character table and the dimensions of the 
conjugacy classes, respectively. Finally, lines 6-12 give the explicit form of the matrices 
for all elements and for all representations of the group. 

The first column of the character table gives the dimensions of the representations. We 
follow the common practice of denoting the representations by their dimensions and 
using primes or numbers to distinguish different representations of the same dimension: 

1, 1', 1", l'", 1 (4) , 1 (5) , 1 (6) , 1 (7) , 1 (8) , 3, 3', 3" (2) 
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Field 


SU(2) L X U(l)y 


u(i)* 


^4 


c 3 


A 4 x C 3 


L 


(2,-1) 


1 


3 


UJ 


3' 


e 


(1, 2) 


1 


1 


u 2 


V 




(1, 2) 


1 


1" 


u 2 


1(8) 


T 


(1, 2) 


1 


V 


u 2 


1(5) 


K 


(2, 1) 





1 


1 


1 


h d 


(2,-1) 





1 


1 


1 


ip T 


(i, o) 





3 


1 


3 


<-Ps 


(i, o) 





3 


UJ 


3' 




(i, o) 





1 


UJ 


1" 



Table 2: Particle content and charges for the model given in ref. [16]. The last column gives the family 
symmetry charges in our notation. In ref. [16], there is an evident typo in the charge assignments to fi 
and r in Section 4 as compared to Section 3 in the same publication. 

Note that we deviate from the notation of ref. [16] where the transformation properties 
of the representation under the factor subgroups are indicated, e.g. 3 <g> u, where co is 
the primitive third root of unity. The reason why we choose another notation is that we 
would like to deal with all groups on equal footing. It is easy to establish the connection 
between the two notations by comparing the representation matrices of A 4 and A4 x C3 
e.g. for the 3: The first, third and fourth generator of A 4 x C3 are identical to the three 
generators 1 of A4 and the second generator generates C3. We can now easily identify 
3 ~ 3 <g) 1, 3' ~ 3 <g> uj, 3" ~ 3®o; 2 . The other cases are handled in a completely 
analogous way (see Tab. 2 for the complete list). Strictly speaking, though, making this 
connection is not necessary. 

The particle content of the model is given in Tab. 2. In the following we list the terms 
that (i) are invariant under the Standard Model gauge symmetry, the i?-symmetry and 
the family symmetry, (ii) contain exactly 2 leptons, (iii) have mass dimension smaller 
than or equal to 6, (iv) are at most linear in the flavon vevs: 

L Lh u h u (p s + L Lh u h u £ + Leh d ip T + L fih d (p T + Lt h d ip T (3) 

To check invariance under the family symmetry we need the decomposition of tensor 
products into irreducible representations (see e.g. ref. [28]) that is readily obtained from 
the character table and the dimensions of the conjugacy classes. E.g. for the first term 
in Eq. (3) we have: 

3 / <g>3 / <g>l®l<g>3 / = (l' + 1 (5) + 1 (8) + 2 x 3")®3' = 2xl + 2xl" / +2xl( 4 ) + 7x3 (4) 

The tensor product contains 2 singlets and thus there are 2 ways to contract the family 
indices to obtain invariant combinations. To do this, however, we need to know the 

^dn the presentation we have chosen, A4 is given by three generators. One can get more information 
on A4 by running the GAP script on the preceding page with the GAP ID [12,3]. 
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Clebsch-Gordan coefficients for A4 x C3, and to our surprise, the general method for 
the calculation of Clebsch-Gordan coefficients for any finite symmetry group is not well- 
known. That is why we have dedicated Appendix D to discuss an algorithm [24] for the 
calculation of Clebsch-Gordan coefficients for finite groups. The first term in Eq. (3) 
after contracting the family indices becomes: 

-^=L 2 L 3 h u h u Lf S ,i + ^7^ 3 Ll hu hu ^,2 + ~^ Ll L z hu h u <fs,3 + ~f^ Ll Ll hu hu € 

+ -^=L 2 L 2 h u h u P, + -^=L 3 L 3 h u h u P, (5) 



After contracting the SU(2) indices and substituting the vevs 

(W) = (hf) = 0, (h^)=v u , (h^) = v d , M = (vt,v t ,v t ), ( W > = («5,0,0) > (0 
Eq. (5) reads: 



(6) 



1 T W T (1) v v + 1 T {1) T {1) v 7; P + 1 T {1) T {1) v v P + 1 T {1) T {1) v 11 P (7) 
— J^ Lj 2 ^3 v u v uVs -r -^-^1 L>\ V u V u ^ + —j=1j2 L, 2 V u V u £ "h — j= -^3 ^3 V u V u t; {( ) 

Following these steps for all the terms in Eq. (3) yields the mass matrices for the charged 
leptons and the neutrinos: 



4 2 7 



L 



(2) 
J 2 

(2) 



e 
1 

v/3 
1 

v/3 
1 



1 

~V3 
1 

2V3 
1 

2^/3 



T 

1 

2V3 
1 

2V3 / 







r(l) 
^2 


4" 




/ J_ 





\ 


If) 





1 


1 


v/3 


2^ 


4 X) 


I 


1 

2^ 


7s / 



The singular value decomposition (here in the special case where the number of rows is 
equal to the number of columns) 



M e+ = D L M e+ D l R , M v = U L M u Uk, 



(9) 



allows us to express the mass matrices as a product of a unitary matrix, a diagonal 
matrix with non-negative real numbers on the diagonal, and another unitary matrix 
where 



D, 



-0.5774+i 0.0000 -0.5774+i 0.0000 -0.5774+i 0.0000 

0.5738-i 0.0636 -0.2319-M 0.5287 -0.3420-i 0.4652 
0.5731-i 0.0702 -0.3474-i 0.4612 -0.2257+i 0.5314 



0.0000 -0.7071 -0.7071 
-1.0000 0.0000 0.0000 
0.0000 0.7071 -0.7071 



(10) 



The neutrino mixing matrix is by definition 



PMNS 



0.8165 + i 0.0000 0.5774 + i 0.0000 0.0000 + i 0.0000\ 
0.4058 -i 0.0449 -0.5738 + i 0.0636 0.0778 + % 0.7028 
0.4052 -i 0.0497 -0.5731 + i 0.0702 -0.0860 - i 0.7019/ 



'IV. 
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In Section 4.8, we discuss how to extract the mixing angles and phases from the most 
general form for £/pmns- in the present case, we obtain 

12 = 0.00, e 2 3 = 35.26, #i3 = 45.00, 5 = -90.00, (12) 

which is tribimaximal mixing. 

Several remarks are in order, (i) At first glance, t/pMNS i n Eq. (11) bears little resem- 
blance to the Harrison-Perkins-Scott matrix {/hps- Note, though, that rephasing the 
fields, Eq. (11) can easily be brought to £/hps form. In other cases, there may also be 
ordering ambiguities (see Section 4.8 for more details), (ii) In contrast to ref. [16] where 
the matrices for the generators in the 3-dimensional representation were wisely chosen 
so that Mi+ is diagonal, our choice for the generators leads to a non-diagonal charged 
lepton mass matrix, see Eq. (8). We have checked that after changing to a basis where 
Mi+ is diagonal (which corresponds to a redefinition of the charged lepton fields), our 
expressions for the Lagrangian, the mass matrices and C/pmns coincide with those in 
ref. [16]. (iii) The same change of basis maps our vevs in Eq. (6) to those of ref. [16]. 
(iv) The reason why our intermediate results do not coincide with those in ref. [16] is 
that we started out with a different choice of generators. Our generators 7\, T 2 , T 3 of 
Ai are connected to the latter ones by S H- T\T 2 T\ and T h> T 2 . (v) For later refer- 
ence, we summarize in Tab. 3 the symmetry breaking patterns for the model at hand. 
In Appendix C on page 32 we will discuss in detail how to find all inequivalent vevs that 
break to different subgroups of a given symmetry. 



1 VEV 


2 VEVS 


3 VEVS 




1)22 — C*3 x C3 


<Pt, £ 


f)io — C3 


<Pt, <fs, € 


fh = l 




1)4 =C 2 




f)4 =C 2 




f)28 = A 4 




1)1 =1 



Table 3: The subgroups to which A4 x C3 is broken when the flavon fields ipx, <ps, £ transforming as 3, 
3', 1", respectively, acquire the vevs (tpx) = {vt,vt,vt), (<fs) = ( u s>0, 0), (£) = v%. The numbering 
of the subgroups corresponds to the output of the GAP script on p. 32. 

Note that since we will extend this analysis to 76 groups, we have to work with the 
generators that are supplied by GAP. It is not feasible to look for the optimal set of 
generators for each of the 76 groups that we will consider; in any case, the results 
are the same, and that the intermediate expressions may be more complicated is not 
relevant, since we have automated the calculation. 

We have written Python programs that interact with GAP to get the generators, the 
character table, the dimensions of the conjugacy classes and the explicit form of the 
matrices for all representations. From this, our code builds the Lagrangian that is 
invariant under all the symmetries, breaks the family symmetry, collects the terms that 
contribute to the charged lepton and neutrino mass matrices, and finally calculates the 
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mixing matrix, the mixing angles and the phases. In Section 4, we will explain the 
details of our systematic scan. 

4. Systematic Construction of the Models 

In the following, we will consider the most general lepton sector with Standard Model 
particle content and up to three flavon fields. For clarity, we summarize our approach 
in form of a flow chart in Fig. 1 on the next page and elaborate on the details in Section 
4.1-Section 4.8. 

4.1. Particle Content 

The particles and their Standard Model charges are listed in the first two columns 
of Tab. 2 on page 6 and will not be reproduced here. To avoid any misunderstandings, 
we emphasize that from now on £ is on the same footing as ipT and <fis and that its 
naming is simply a relic from earlier sections. 

We restrict ourselves to such models where the lepton doublet L transforms in a 3- 
dimensional representation and e, fi, r transform in 1-dimensional representations of the 
family symmetry. Plausible as this may sound, there is no physics reason for that, but 
rather, as we will explain below, without these assumptions the number of family charge 
assignments quickly grows too large to allow for a systematic scan. 
Regarding the Higgs sector, we will assume that there are exactly two fields. For one 
thing, we have supersymmetric models in mind that require an even number of Hig- 
gses. For another, more than two Higgs fields would spoil the unification of the gauge 
couplings. Thus, h u , hd are assigned any 1-dimensional representation. 

4.2. Family Symmetry 

In the list of all groups of order < 100 [23], we find 90 groups which have a 3-dimensional 
representation (see Tab. 4 on page 24). We iterate over 76 out of these 90 groups that 
can be scanned in less than 60 days. Note that abelian groups only have 1-dimensional 
representations, and are thus not included in our scan. For a systematic scan of discrete 
abelian symmetries, see ref. [29]. At this point in the algorithm, we calculate the relevant 
information on the group that we will need later on: (i) The dimension of the group, 
the number and dimensions of its conjugacy classes, and its character table, (ii) the 
irreducible representations and their tensor products, (iii) the representation matrices 
for the irreducible representations, (iv) the Clebsch-Gordan coefficients depending on 
the choice of the representation matrices. 

4.3. Charge Assignments 

We iterate over the inequivalent family charge assignments. As mentioned before, L 
is assigned any 3-dimensional representation, and e, /i, r, h u , hd are assigned any 1- 
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Choose group G; cal- 
culate CG-coefficients 







No more charges 



Assign discrete charges to 

L, e, fx, T, ft*, k d , <fT, ^S: £ 







Find invariant terms with 2 leptons; 
couplings LL, Le, Lfi, Lr present? 




yes 


Contract discrete symmetry in- 
dices; Lagrangian already known? 




new Lagrangian 


Remember Lagrangian; 
contract SU(2) indices 



Next vev 



Repeat choosing vevs; Charged 
lepton mass matrix singular? 


No vevs left 






no 


Calculate neutrino mass ma- 




trix; calculate Ufmks ail d 




Euler angles; write to file 





Figure 1: Systematic scan for models with family symmetry G and up to three fiavon fields. 



dimensional representations. We do not make any assumptions on the representations 
of the fiavon fields <pr, <£s, £• 

In absence of a mechanism for generating mass hierarchies, we cannot distinguish be- 
tween e, /i and r; they have the same quantum numbers and their naming is largely a 
matter of convention. Thus, to avoid iterating over configurations that give the same 
physics, we consider any permutation of the charge assignments to the e, /i, r to be 
equivalent. Once the mixing matrix has been derived, we can reorder its rows to recover 
the cases corresponding to the aforementioned permutations. In other words, we identify 
the electron, muon and tau a posteriori and rename them where necessary. 
The same holds for the fiavon fields ipx, <Ps, £5 an d since their interactions are not directly 
observable, it is not even necessary to rename them. 

The running time for the algorithm scales with the number of irreducible representations. 
Let us denote by N\, N3 and N a the number of 1-dimensional, 3-dimensional and all 
irreducible representations, respectively. Then the total number of inequivalent family 
charge assignments is 

N 3 x £.(#1,3) x iWiVx x C r (N a ,3), (13) 

where the first factor corresponds to L, the second one to e, /1, r, the third one to h u , hd, 
and the last one to <px, <Ps, £• Consider e.g. the second factor. The charge assignments 
to e, /1, r do not depend on their order, so if all three charges are distinct, the number 
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of inequivalent choices is given by iVi-choose-3. In the more general case, two or more 
charges may be the same, and the number of inequivalent choices is given by (for N = N\ 
and k — 3): 

(AT + h - l)\ 

C r (N, k)=C(N + k- 1, k) = l (iV _ 1)!A | (14) 

The group A4 x C3 has N a = 12 irreducible representations with A3 = 3 and N± = 9. 
Thus, the total number of inequivalent charge assignments is 14,594,580. 

Note that the case of less than three flavons is automatically included in the algorithm, 
since for a given flavon field we are also iterating over the vev v = which effectively 
removes the corresponding field from the Lagrangian. 

4.4. Invariant Lagrangian 

In principle we could now construct the most general Lagrangian that is invariant under 
the gauge and family symmetries and contract the family indices. As a matter of fact, 
that is what we had initially done. However, because of the large number of inequivalent 
charge assignments, it is a better approach to try to determine as early in the algorithm 
as possible whether a given model is viable or not. 

For deriving the mass matrices of the charged and neutral leptons, we only need those 
terms in the Lagrangian that contain exactly 2 leptons, at most 2 Higgses and at most 
1 flavon, since terms that have mass dimension greater than 6 or that are quadratic 
in the flavon fields are suppressed. We establish the invariance of a given term under 
the (gauge or family) symmetry by checking whether the tensor product of the particle 
representations contains a singlet. Note that this operation is very "cheap" for the 
computer as compared to doing the full contractions. 

If there are no invariant terms at all or some of the couplings that we need for giving 
masses to the leptons are absent, we can skip the rest of the calculation and immediately 
continue with the next assignment of family charges. It is important to note that such 
improvements to the algorithm are crucial for keeping the running time within reasonable 
limits. 

4.5. Contracting the Indices 

To contract the family indices, we use the Clebsch-Gordan coefficients that we have 
already calculated in the first step of the algorithm. Their derivation for an arbitrary 
finite group (and choice of representation matrices) is not well-known, and the results 
available in the literature cover only specific cases. 

One algorithm for the general case that we are aware of was presented in ref. [21]. We 
have implemented an algorithm due to van den Broek and Cornwell [24] that we believe 
to be more efficient and that is discussed Appendix D. 

After all gauge and family indices have been contracted, we rearrange the terms in the 
Lagrangian and the particles in each term to bring them into lexicographical order. 
We then compare the Lagrangian at hand with the list of Lagrangians from previous 
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iterations. If the given Lagrangian is already known, we continue with the next iteration 
over family charges. Otherwise, we save it to the list and contract the SU(2) indices. 
For A4 x C3, the number of different Lagrangians is 39,900. 

4.6. Substituting the Vacuum Expectation Values 

In Appendix C we show how to find all vevs that break to a particular subgroup of the 
flavor symmetry. Unfortunately, we have examples which show that two different sets 
of vevs that induce the same symmetry breaking pattern can lead to different values for 
the mixing angles, so our classification of the vevs does not help us in classifying the 
models. 

In lack of a better approach, we choose the entries of the vevs to be or 1. For A 4 x C3, 
the number of different vev configurations per family charge assignment may range from 
2 3 to 2 9 (<pt, tpsi £ transform all in a 1-dimensional or 3-dimensional representation, re- 
spectively). Then we replace all fields by their vacuum expectation values. Analogously, 
we substitute hff = v u , hp = Vd and hu = h d = 0, where the superscripts denote the 
SU(2) indices. 

4.7. Mass and Mixing Matrices 

By construction, the Lagrangian contains only terms quadratic in the lepton fields whose 
coefficients give the charged and neutral lepton mass matrices and M v , see Eq. (8) 
on page 7. The singular value decomposition diagonalizes the mass matrices by unitary 
transformations whose product give the neutrino mixing matrix: 

M e+ =D L M e+ rf R , M u = U L M u Ul Upmns = D L Ul (15) 

If the charged or neutral lepton mass matrix is singular (i.e. at least one of the masses 
is zero), we continue with the next iteration over the vevs. Only when all of the vevs 
are exhausted, we continue with the next iteration over the charge assignments. 

4.8. Euler Angles and Recognizing Tribimaximal Mixing 

Of the many different parametrizations [30] for the neutrino mixing matrix that are 
mathematically equivalent and describe the same physics, we follow the standard nota- 
tion as advocated by the Particle Data Group [31]. To extract the mixing angles and 
phases, we use the explicit formulae presented in ref. [27] that use a slightly different 
convention for the two Majorana phases. 

Note that t/pMNS = D^Ul is guaranteed to be unitary by virtue of the singular value 
decomposition, but (before using the rephasing freedom) may not necessarily be in the 
standard form as given in ref. [31]. Luckily, the formulae in ref. [27] are applicable for 
any unitary matrix so that we can circumvent this technical complication. 

As explained before, for a given charge assignment to e, fi, r, we do not iterate over all 
its permutations, since we have the freedom of renaming the particles. As a consequence, 
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when calculating the mixing angles, we must consider all permutations of the rows of 
U-pmns- The naming of the neutrinos is then fixed by the corresponding charged leptons. 

4.9. Details of the Technical Implementation 

The programs were written in Python 2.6.2 [32] and used Numpy 1.4.1 [33] for the 
linear algebra operations. For the group theory calculations, we interfaced our programs 
with GAP 4.4.12 [20]. For data analysis we used PyROOT [34] that provides a convenient 
interface to Root 5 . 17 [35]. For generating the graphs, we mainly relied on rootplot [36], 
but in some cases we had to extend its functionality by overloading its classes and directly 
using MatPlotLib 1.0.0 [37]. The code was executed on the Linux cluster at the Centre 
de Calcul de I'Institut National de Physique Nucleaire et Physique des Particules in 
Lyon, France. 

After the first initialization (creating the data on the group, opening files for reading 
and writing, etc.), it takes less than 1 second to calculate the full details for the model 
corresponding to the irrep assignments in Section 3 and 128 choices of vevs. As expected, 
we rediscover tribimaximal mixing (for six vev configurations). 

4.10. Running Times for the Different Groups 

For Ai x C3, the execution time was 16 hours and 43 minutes on one computer with a 
3 GHz Intel Xeon processor. From this, we get a useful measure to assess the running 
times for the other groups in Tab. 4 on page 24, since we know the number of inequivalent 
charge assignments for each group before we start running the programs. Note, however, 
that this can give us only an order of magnitude, because the computer center does not 
guarantee the same hardware on all of its machines. 

The running times for the 90 groups of order < 100 that have a 3-dimensional represen- 
tation ranges from 29 seconds to ca. 177 years, whereas they do not necessarily increase 
with the number of elements in the group, but depend on the number of representations. 
For the present publication, we have decided to limit ourselves to those 76 groups that 
can be scanned in less than 60 days on one computer and have indicated the 14 groups 
that exceed our time limit by red text color in Tab. 4. For the actual calculations, we 
have distributed the calculations on more than one computer. 

It should be noted, though, that it is not impossible to tackle the groups that require 
longer running times. For one thing, one can add more theoretical priors to reduce the 
number of configurations that need to be scanned. For another, one can rewrite time- 
critical components in C++ and integrate them into the Python programs. Yet another 
way to improve the code is to use parallel computing. 

5. Phenomenology 

We will now discuss the results that we obtained from the systematic scan of family 
symmetries, charge assignments and vacuum configurations. We refer the reader also 
to Tab. 4 on page 24 where some of the results of this section will be summarized. 
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It is important to stress that we are not specifically searching for tribimaximal mixing, 
but constructing all models for a given symmetry group (with the qualifications detailed 
in Section 4.3 and Section 4.6). 

We will only list inequivalent models: We consider two models to be equivalent, if 
their Lagrangians after contracting the family indices, but before symmetry breaking 
are equal. In the plots, however, the data points correspond to vacua and some may 
correspond to the same Lagrangian. 

Since we will be discussing many different groups that may not all have a standard name, 
we will use their GAPIDs, and e.g. denote A 4 x C 3 by 0(36,11). The correspondence 
between the GAPIDs and the groups is given in Tab. 4 on page 24. 

5.1. Models for A 4 x C 3 

We will start with the results for the "classic" group 0(36,11) = A A x C 3 . The 
14,594,580 family charge assignments to L,e, /i,T,h u ,h ( i,(pT, l Ps,^ give 39,900 inequiv- 
alent Lagrangians out of which 22,932 have non-singular charged lepton and neutrino 
mass matrices. 

In this set, we find 4,233 models of tribimaximal mixing (18.5%). For a given model, 
there may exist more than one vacuum configuration that leads to TBM (e.g. 6 vevs for 
the model in Section 3), and we have not counted them separately. 
4,481 models (19.5%) lie in the 3a range of their measured values. We find no models 
that lie in the la range, because 9 i3 = 0° is excluded at la (cf. Tab. 1 on page 5). This 
fact is nicely illustrated in the third plot of Fig. 2(b): The la range, represented by the 
green band, is empty. 

Fig. 2 shows the distribution of the mixing angles #12, #23, and #13, where we are now 
counting the vacua and not the models. The reason for this is that for one and the 
same Lagrangian, the values of the mixing angles will in general depend on the choice 
of vacuum. The histograms in Fig. 2(a) have 15,992,118 entries reflecting the fact that 
for each of the 39,900 inequivalent Lagrangians, we are looping over 8 to 512 vacua, 
depending on the dimensions of the irreps assigned to the flavon fields. 
The area of each histogram has been normalized to 1 and the bin width is 1°, so the y- 
axis gives the percentage of vacua that produce the angles on the x-axis. The green and 
yellow bands correspond to the la and 3a ranges, respectively (cf. Tab. 1 on page 5). In 
Fig. 2(a), we simply count the number of times that a given angle is reproduced irrespec- 
tive of the values that the two other angles may take. E.g. from the first histogram we 
can read off that 7.6% of the vacua give a value for dyi that is consistent with experiment 
at 3o", where #23 and 6*13 can take any values. 

Now we investigate whether we can obtain some predictions by introducing priors. In 
Fig. 2(b) we have restricted 2 of the angles to their 3a intervals and plotted the third 
one. As a consequence, the numbers of entries in the histograms are not equal. 
The most striking difference between Fig. 2(a) and Fig. 2(b) is that now the most likely 
value for 6*23 is 45° and lies in the la interval. Furthermore, the number of vacua in the 
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(a) Number of models that give 9ij with no constraints on the other 2 angles. Each histogram has 
15992118 entries. 
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(b) Number of models that give Oij with the other 2 angles restricted to their 3er interval. The histograms 
have 838289, 148886 and 225844 entries, respectively. 



Figure 2: Number of models with family symmetry ©(36, 11) = A4 x C 3 that give the mixing angles 
denoted on the x-axis. The area of the histograms is normalized to 1 and the bin width is 1. The green 
and yellow bands correspond to the la and 3ct ranges, respectively. 



experimentally disfavored region has decreased significantly, and 58% of the vacua are 
in the 3a interval. 

For #13, values near 90° are now excluded, and the 3a interval is almost depopulated 
except for #13 = 0°, which at the same time corresponds to the maximum of the his- 
togram. This can be interpreted as a prediction for #13 to be 0° (at leading order) based 
on current experimental data and the theoretical assumption of an A± family symmetry. 
38% of the vacua are in the 3a interval. 

For #12, the most likely value is still 0°, but 35° is now the third- most assumed angle. 
Clearly, the experimental data on #23 and #13 is pushing us in the right direction. 11% 
of the vacua are in the 3a interval. 

It is an interesting observation that the effect of our "cuts" were such that the preference 
of the data for the experimentally allowed ranges became much more pronounced. This 
is an unexpected and non-trivial result and may further testify to the phenomenological 
viability of A±. 

To learn more about the correlation of the angles and how priors may affect them, in 
Fig. 3 we present the distribution of 2 out of 3 angles, respectively. The color bar on the 
right-hand side of each figure gives the correspondence between the colors in the plots 
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(b) Number of models that give Oij and 9 mn with the remaining angle restricted to its 3<r interval. The 
histograms have 2941000, 3675600 and 1057170 entries, respectively. 

Figure 3: Logarithmic plot of the number of models with family symmetry ©(36, 11) = A4 x C3 that 
give the mixing angles denoted on the axes. The bin width on both axes is 1. The base of the logarithm 
is 10, and the color map on the right side of each plot gives the exponents. 

and the logarithm to base 10 of the number of vacua with the angles 0y and 9 mn on the 
x- and y-ax.es, respectively. 

In Fig. 3(a), each 2-dimensional histogram has 15,992,118 entries which correspond to 
the full set of vacua that we had also previously considered in Fig. 2(a). In analogy 
to Fig. 2(a), we have imposed no constraints on the third angle that is not plotted. 
From the first plot, we cannot read off much, except that there exist certain "hot spots" 
(e.g. 9 12 = 0° and 6*23 = 45°) that correspond to large numbers of vacua, and that the 
regions near the lower corners are by comparison less populated. In the second and 
third plots, we see that there are considerably fewer models for 6*12 > 70° and 6*23 > 75°, 
respectively. In the case of the second plot, this holds even for much lower values of 
9±2 ^ 35°, given that #13 is not larger than ~ 15° or near 0°. 

In Fig. 3(b), we present the same correlations as in Fig. 3(a), but this time, we have 
required that the third angle be in its 3<r interval. As a consequence, the numbers of 
entries in the histograms are not equal. We have used the same normalization of the 
color bars in Fig. 3(a) and Fig. 3(b) to facilitate comparisons between them. 
Considering the first plot, we can see that most vacua lie in a band 6*12 — 30° — 60°, 
whereas this effect becomes less pronounced for #23 — 15° — 30°. For #23 — 45°, values 
of 9 12 = 0°, 35°, 55°, 90° are favored (red bins in plot). 
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The second plot clearly shows that 9 12 > 70° and to a lesser extent 9 12 < 10° are 
disfavored. For 9i2 — 10° — 70°, a band of #13 = 0° — 10° that widens with increasing 
9 12 is sparsely populated, but note that for 9 Vi = 0° and 9 12 = 45°, there is one of the 
highest counts of vacua in the plot as indicated by the red bins. This is consistent with 
our previous observation from Fig. 2(b) that 6* 13 = 0° is preferred, but that otherwise 
the region #13 < 10° is disfavored. 

In the third plot we again observe a band structure #23 — 15° — 60° where most of the 
vacua are concentrated, and find that #13 < 10° and #13 > 80° are disfavored. The 
combination #23 ~ 45° and #13 = 0°, however, is not preferred. As indicated by the red 
bins, the most likely combination of angles lies elsewhere. 



Figure 4: Scatter plot for #12, #23 ? $13- Each point corresponds to a bin in the 3-dimensional histogram 
that has at least 1 entry. The bin width is 1. The color of the points (from blue to red) correspond to 
the logarithm of the number of models (from lower to higher) . 

Ideally, to give a graphical representation of the full information on the angles and their 
correlations, we would create a 3-dimensional histogram with l2 , #23, and #13 on the 
axes, and plot the number of vacua along a 4th dimension. Since this is not feasible, we 
present a plot in 3 dimensions, where the color of the data points indicates the number 
of vacua. 

In Fig. 4(a), each point represents a bin in a 3-dimensional histogram: If there is at 
least one vacuum that produces the angles (6*12, 6*23, 6*13), we set a point at the respective 
coordinates. The bin width on each axis is 1, and in total there are 90 x 90 x 90 bins, 
of which 5,528 are not empty. The color of the point denotes the logarithm to base 
10 of the number of vacua that give the respective angles, where the colors from blue 
to red correspond to an increasing number of vacua. We have not displayed the color 
map for the plots, since we find it difficult to extract quantitative information from the 
3-dimensional representation and rather use it as a means of uncovering correlations 
between the angles and the qualitative features of A 4 as a symmetry group. 
In Fig. 4(b), we display only those 1,287 bins that have more than 1,000 entries. This 
removes much of the cluttering and gives a clearer picture regarding where the most 
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(a) The 5,528 bins that are > 1. 



(b) The 1,287 bins that are > 1000. 
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likely vacua are concentrated. One immediate observation is that the parameter space 
for (#i2, 023, #13) is not uniformly populated: There are very few vacua for 612 — 60° — 90°; 
low and high values of #13 are disfavored (except 6>i 3 = 0°); and rotating Fig. 4(b) to 
view it from different perspectives, we see that most of the vacua are concentrated in a 
volume 12 ~ 10° - 60°, 23 ~ 20° - 70°, 13 > 15°. 

Considering the $12 — 023 plane of Fig. 4(b) that corresponds to 6*13 = 0°, we find that 
012 = 023 = 45° are the most likely values (see red points), which is in agreement with 
the first plot in Fig. 3(a). 

Regrettably, we fail to see any preference for tribimaximal mixing or the experimentally 
allowed values. Without putting in at least some priors, the best we can do is setting 
approximate upper and lower bounds on the angles. 




g 0.04 - 




180 



(a) Number of models that give the phase denoted on the cc-axis with no constraints on 612, 623, 6*13. 
Each histogram has 15992118 entries. 




S 0.10- 




(b) Number of models that give the phase denoted on the x-axis with 612, O23, #13 restricted to their 3a 
intervals. Each histogram has 86014 entries. 



Figure 5: Number of models with family symmetry (5(36, 11) = A 4 x C3 that give the Dirac and Majo- 
rana phases denoted on the x-axis. The area of the histograms is normalized to 1 and the bin width is 
4, so one must multiply the height of the bars by 4 to get the percentage of vacua. 



If 0i 3 7^ 0°, the CP-phase S becomes relevant. In Fig. 5, we present the distribution 
of the vacua for 5 and the two Majorana phases <p\ and $2 (for our conventions, see 
ref. [27]). Since the phases may take values from -180 to +180, we have chosen a bin 
width of 4. The area of the histograms is still normalized to 1, but now one has to 
multiply the height of the bars by a factor of 4 to get the percentage of vacua. 
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In Fig. 5(a), we simply count the number of vacua that realize the phase denoted on the 
x-axis, irrespective of the values that the other 2 phases and the 3 mixing angles may 
take. In Fig. 5(b), however, we use the experimental information that is available to us, 
namely, we restrict 6*12, #23 and 6^3 to their respective 3a intervals; experimental data 
on the phases is not available. 

We get a clear prediction that 5 = 0i — 02 = 0, which changes to 5 = 02 = and 
0i = ±90 if we take the experimental constraints into account. One reason for this 
may be that for the flavon fields, we have chosen only real vacuum expectation values. 
Complex numbers are introduced in the Lagrangian only through the Clebsch-Gordan 
coefficients. 

In analogy to Fig. 3 and Fig. 4, we have also analyzed the 2-dimensional histograms and 
the 3-dimensional scatter plots, but abstain from reproducing the graphs in the present 
publication. From the 2-dimensional histograms with the experimental constraints im- 
posed, we learn that the phases prefer a discrete sets of values, and one may get hints 
at some correlations, e.g. 02 = 0i ± 180. The 3-dimensional scatter plot is consistent 
with a random, uniform distribution of the phases, but we can distinguish some lines 
corresponding to a higher concentration of vacua. 

5.2. A Model with 6 13 ~ 5 

Recently some analyses have reported on possible hints of a non-zero #13 [6, 17, 18]. 
A model with 0° ^ 9 13 < 8.6° that lies in the la interval is easily constructed. For 
the family symmetry, we take 25(36, 11) = A4 x C3 and assign the family charges 
(L, e, fj,, r, h U7 h d , (p T , <p s , f) ~ (3, 1, 1 (4) , 1'", 1', 1 (7) , 1 (5) , 3", 3"). When the flavon fields 
acquire vevs along the directions ((pr) = (1), ((fs) = (1,0, 1), (£) = (1, 1, 1), we obtain 
the mixing angles 9 12 = 33.9°, #23 = 40.9° and 6*13 = 5.1° that all lie in the la interval of 
the experimentally determined values. Incidentally, we have chosen the model in such 
a way as to produce a 9 ri that is close to the present best-fit value of ref. [6] with the 
modified Gallium cross-section. 

5.3. Results for 76 Flavor Groups 

In this section, we present an overview of the models that we obtained from our system- 
atic scan. Due to the sheer volume of the data, we will limit ourselves to making a few 
qualitative observations and relay the detailed analysis to an upcoming publication. 
One main result of our analysis is that we have found thousands of new models that 
give exact tribimaximal mixing. Fig. 6 shows the number of inequivalent models for 
each of the 76 groups that (i) are of order < 100, (ii) have a 3-dimensional irreducible 
representation, and (iii) can be scanned in < 60 days (see Section 4.10 for more details). 
We have excluded those Lagrangians from our analysis that lead to a singular neutrino 
or charged lepton mass matrix. The red bars indicate the fraction of Lagrangians for 
which at least one choice of vevs leads to tribimaximal mixing, and the blue bars give 
the number of models that lie in the 3a interval of the measured angles (cf. Section 2 
on page 4). The correspondence between the GAPIDs on the x-axis and the full name of 
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Figure 6: The number of models per symmetry group. On the x-axis, we label the flavor symmetry 
g by its GAPID, cf. Tab. 4 on page 24. The red and blue color of the labels on the x-axis indicates 
that q D A4 and g C U(3), respectively, whereas the green color signifies that both conditions are 
satisfied simultaneously. Along the y-axis, the blue bars give the number of Lagrangians that lead to 
non-singular mass matrices. The green bars indicates the number of models that lie withing the 3cr 
interval, and the red bars finally give the number of models for which at least one vacuum configuration 
gives tribimaximal mixing. 



the group is given in Tab. 4 on page 24, where we also list the exact numbers of models 
that may be difficult to read off from the graph. 

The conspicuous gaps in Fig. 6 are a consequence of our criterion that the mass matrices 
be non-singular, i.e. we do not consider such cases where any of the neutrinos (or charged 
leptons) is massless. 

Out of the 76 groups that we scanned, 9 (12%) have only singular mass matrices. 44 
groups (58%) lie in the 3a interval, and 38 (50%) are even tribimaximal (for at least one 
vacuum configuration, respectively). Only for 23 groups (30%) we could not find any 
vacuum configuration that satisfy the experimental limits. Note, though, that despite 
being very general, our scan is not fully comprehensive, since (i) we assume that the 
lepton doublet transforms in a triplet, and (ii) we do not scan over all possible vevs. 
Owing to this fact, there may exist even more viable models than we have identified. 
To explore the connection between tribimaximal mixing and A4, we have color-coded 
the group names on the x-axis of Fig. 6. The blue, red and green colors correspond to 
C U(3), D A*, and A 4 C fl C U(3), respectively. Of the 76 groups, 35 contain A 4 as 
a subgroup, but only for 16 out of these 35 groups we can find vacua that give models 
of tribimaximal mixing. It is conceivable, though, that one may find TBM models for 
the other groups, if one introduces more than 3 flavon fields. 

The chances of finding TBM does not scale with the total number of models, as e.g. (5(81, 7) 
and the next adjacent three groups show. Yet again, this conclusion may heavily depend 
on the number of flavon fields. 

We find four groups, 0(84,9), (5(96,198), (5(96,201) and (5(96,202) that are neither 
subsets of U(3) nor contain an subgroup and nevertheless can accommodate models 
of tribimaximal mixing. 

An interesting observation from Fig. 6 is that for 9 groups the models with TBM are 
identical to those in the 3a interval (the green bars are almost completely covered by 
the red ones). 5 of the 9 groups belong to the T-series of SU(3) subgroups [38]. 
The smallest group for which we find TBM is (5(21, 1) = T7 which after A4 is the smallest 
group that has a 3-dimensional irreducible representation. The family symmetry T7 has 
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Figure 7: Number of models for ©(21, 1) = T7 that give 9ij with the other 2 angles restricted to their 
3<T interval. The histograms have 4472, 984 and 1920 entries, respectively. The area of the histograms 
is normalized to 1 and the bin width is 1. The green and yellow bands correspond to the la and 3tr 
ranges, respectively. 



been studied in refs. [39-41]. 

It is also worth mentioning that the second smallest group with models of TBM is 
0(24, 13) = Ai x C2 where the C3 factor of the model in ref. [16] has been replaced by 
a C 2 . 

In Fig. 7 we show for 0(21, 1) = the distribution of the mixing angles. If we use the 
experimental data on 2 of the angles and plot the multiplicity of the third one, we find 
that 6*12 — 35° is the second-most likely angle to be produced, and the only one within 
the 3cr interval. For 623, we obtain a unique prediction 823 = 45°. As for #13, the value 
0° is both the most likely angle as well as the only one attained within the 3a interval. 

As compared to 0(36, 11) = A^x C3, there are other groups that have a larger fraction of 
TBM models. Consider e.g. 0(39, 1) = C13 C3 = T13, where we find 288 inequivalent 
models, of which 171 (59%) are TBM. Remarkably, any model for that group that has a 
vacuum for which the mixing angles are consistent with experiment at the 3a level also 
allows for TBM to be realized. 



6. Conclusions 

In this publication we scanned 76 groups and constructed a total of 439,820 Lagrangians 
out of which 59,019 are consistent with experiment and 31,137 are tribimaximal. The 
large set of viable models allowed us to look at correlations between the mixing angles 
and make a prediction for # 13 that will be measured in upcoming experiments. 

We have presented an explicit model with 13 = 5.1° to show that the recent tentative 
hints of a non-zero #13 can be accommodated. We found tribimaximal mixing in 38 
flavor groups; most of these groups had not been considered for model building before. 
We hope that the calculational tools and methods that we have outlined will be useful 
for future model building efforts. 

We would like to emphasize that we do not advocate a probabilistic approach to model 
building along the lines of the landscape idea in string theory. Rather, we are trying to 
maximize our chances of finding the correct model(s) by starting out with a large set that 
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reproduces the mixing angles within the current experimental limits. In future, we plan 
to take this analysis several steps further and look at the generation of mass hierarchies, 
the vacuum alignment problem, and finally include the quark sector. Invariably, each 
step will drastically reduce the number of models, and the goal is to find at least one 
that passes all criteria. 

On the other hand, for answering the question whether any discrete flavor group is 
inherently connected to tribimaximal mixing, a probabilistic approach may be useful: 
The easier tribimaximal mixing can be realized in a given group, the more pronounced is 
the connection. In this sense, At fares well, but T i3 and maybe T 7 should be considered 
to be on the same footing, if not more promising. 
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A. List of Groups of Order at Most 100 



Only very few groups were given dedicated names by the mathematicians and physicists 
who studied them. Examples are the cyclic groups C n , the symmetric groups S n , the 
alternating groups A n and the dihedral groups D n . The vast majority is described by 
their substructure and a "prescription" of how to put together these parts to form the 
full group. In Appendix E we sketch the concepts that are crucial for the description of 
the groups: the direct product, the semidirect product and the short exact sequence. 

Since the list of 1,048 groups of order < 100 is too long to include in the present 
publication, we have made it available for download [23]. For details on the generation 
of this list, see Appendix B. 

The list of the 90 groups of order < 100 that have a 3-dimensional irreducible repre- 
sentation and that we have systematically scanned for viable models of lepton flavor is 
given in Tab. 4 below. 

A.l. Notation and Conventions 

Our notation follows the GAP Reference Manual p. 356 [42] with the following excep- 
tions. We denote the direct product by "x" and the semidirect by iV K where N is 
normal. Beware that this convention is not unique and that the symbol "x^," may point 
the other way. In writing short exact sequences like 1— > N — > G — > Q — > 1, we will 
omit the leading and trailing trivial groups in order to make our notation more compact. 

We denote the dihedral group of a regular n-gon by D n , and not by D 2n , as some authors 
prefer to do. C n or Z„ is the cyclic group of order n. S n and A n are the symmetric and 
alternating groups, respectively. Q 4 and Q 8 are the quaternion and octonion groups, 
respectively. SL(n,p) is the special linear group over a finite field, i.e. the set of all n x n 
matrices with determinant 1 and values from a field of order p. 

To facilitate the comparison with the existing literature, we list in Section A. 3 a non- 
exhaustive list of alternative names for some of the groups considered in our analysis. 

Many of the groups that we consider do not have specific names, and we will refer to 
them by their GAPIDs. (3(m,n) will denote the group that is generated in GAP by the 
command SmallGroup(m,n). 

A.2. The List of 90 Groups 

In Tab. 4 below, we list the 90 groups of order < 100 that have a 3-dimensional irreducible 
representation. 

The first column gives the GAPID which is a label that uniquely identifies the group in 
GAP. The first number in the square brackets is the order of the group, and the second 
number simply enumerates different groups of the same order. The GAPIDs of the 14 
groups that require more than 60 days of computer time (see Section 4.10) are marked 
in red. 



23 



The second column gives the name of the group. If two or more groups by the same 
name are isomorphic, we list only one. For the conventions we used in naming the groups 
and for a non-exhaustive compilation of alternative names common in the physics and 
mathematics literature see Section A.l and Section A. 3, respectively. 
The third column indicates whether the group G is a subgroup of U(3). If G is in SU(3), 
an orange check mark is shown (i/), otherwise a blue one (i/). 

The fourth and fifth columns indicate whether G is in U(2) or U(2) x U(l), respectively 
(blue check mark). If G is in SU(2) or SU(2) x U(l), respectively, the check mark is 
orange. 

The sixth column indicates whether G contains A 4 as a subgroup. The color of the check 
marks has no significance. 

The seventh column gives the total number of inequivalent models, and the eight and 
ninth columns show the number of models that have vacua with mixing angles that lie 
in the 3a interval or are tribimaximal, respectively. 



Table 4: The 90 groups of order < 100 that have a 3-dimensional irreducible representation. For details, 
refer to the text preceding this table. 



GAPID 


Group 


u 3 


u 2 


U 2x i 


A A 


Models 


3(7 


TBM 


[12,3] 


A 4 


✓ 






✓ 


90 








[21,1] 


CV ~*<p C3 


✓ 








108 


36 


36 


[24,3] 


SL(2,3) 


✓ 








135 








[24, 12] 


Si 


✓ 






✓ 











[24, 13] 


C 2 x A A 


✓ 






✓ 


2034 


344 


288 


[27,3] 


(C 3 x C 3 ) x v C 3 


✓ 








34992 


2430 





[27,4] 


Cg X v C 3 


✓ 








34992 


4536 





[36,3] 


(C 2 x C 2 ) x v C7 9 


✓ 








53535 


10621 


3459 


[36, 11] 


C* 3 x A 4 


✓ 






✓ 


22932 


4481 


4233 


[39, 1] 


C13 X<p C 3 


✓ 








288 


171 


171 


[42, 2] 


C 2 x (C 7 x v C7 3 ) 


✓ 








2682 


445 


273 


[48,3] 


(C 4 x C 4 ) x v C 3 


✓ 






✓ 


270 


90 


90 


[48, 28] 


SL(2,3) -> C 2 


✓ 

















[48, 29] 


GL(2,3) 


✓ 


✓ 


✓ 













[48, 30] 


A4 x v C 4 


✓ 






✓ 


48 








[48,31] 


C 4 x Ai 


✓ 






✓ 


14937 


2864 


2712 


[48, 32] 


C 2 x SL(2,3) 


✓ 








2052 


344 


288 


[48, 33] 


SL(2,3) x v C 2 


✓ 


✓ 


✓ 




2052 


344 


288 


[48, 48] 


C 2 x 5 4 


✓ 






✓ 


16 








[48, 49] 


C 2 x C 2 x A4 








✓ 


5805 


640 


561 


Continued on next page 
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Table 4 - continued from previous page 



GAPID 


Group 


f/ 3 


u 2 


^2x1 




Models 


3(7 


TBM 


[48, 50] 


(C 2 x C 2 x C 2 x C 2 ) x p C 3 








✓ 


189 








[54,8] 


((C 3 x C 3 ) x v C 3 ) x v C 2 


✓ 

















[54, 10] 


C 2 x ((C 3 x C 3 ) x v C 3 ) 


✓ 

















[54, 11] 


C 2 x (C 9 x v C 3 ) 


✓ 

















[57, 1] 


Cig x v C 3 


✓ 








405 


198 


198 


[60, 5] 




✓ 






✓ 











[60, 9] 


C 5 x A 4 


✓ 






✓ 


11575 


2063 


1983 


[63, 1] 


C7 X v Cg 


✓ 








24345 


3792 


795 


[63, 3] 


C 3 x(C 7 x v C 3 ) 


✓ 








15246 


2483 


1863 


[72, 3] 


C/8 X v Cg 


✓ 


✓ 






18714 


3272 


1344 


[72, 15] 


((C 2 x C 2 ) x v C 9 ) x v C 2 



















[72, 16] 


C 2 X ((C 2 X C 2 ) X^Cg) 


✓ 

















[72, 25] 


C* 3 x SL(2,3) 


✓ 


✓ 


✓ 




18105 


3441 


3261 


[72, 42] 


C 3 x S 4 


✓ 






✓ 


108 








[72, 43] 


(C 3 xA 4 )» v C 2 








✓ 











[72, 44] 


Ai x S* 3 








✓ 


2451 


399 


336 


[72, 47] 


C 6 x A 4 


✓ 






✓ 











[75,2] 


(C 5 x C 5 ) x v C 3 


✓ 








477 


234 


234 


[78, 2] 


C 2 x (C 13 x„ C 3 ) 


✓ 








2541 


810 


810 


[81,3] 


(C 9 XC 3 )X V C 3 



















[81,4] 


Cg X^ Cg 



















[81,6] 


C 2 7 X^ C 3 


✓ 

















[81,7] 


(C 3 x C 3 x C 3 ) x v C 3 


✓ 








24329 


1296 





[81,8] 


(C 9 XC 3 )X V C 3 


✓ 








32416 


1782 





[81,9] 


(Cg X C 3 ) X^C 3 


✓ 








32076 


1782 





[81, 10] 


(C 3 X C 3 ) -»(?-> (C 3 X C 3 ) 


✓ 








20736 


1161 





[81, 12] 


C 3 x ((C 3 x C 3 ) x v C 3 ) 



















[81, 13] 


C 3 X(C 9 X ¥ ,C 3 ) 



















[81, 14] 


(Cg X C 3 ) X^C 3 


✓ 

















[84, 2] 


C 4 x (C 7 x v C 3 ) 


✓ 








4752 


714 


567 


[84, 9] 


C 2 x C 2 x (C 7 x v C 3 ) 










2136 


366 


306 


[84, 10] 


C 7 x A 4 


✓ 






✓ 











[84, 11] 


(C X4 x C 2 ) x) v C 3 


✓ 






✓ 


678 


192 


192 


[93, 1] 


C 3 i x,, C 3 










507 


249 


249 


[96, 3] 


((C 4 x C 2 ) x v C 4 ) x v C 3 








✓ 


324 


90 


90 


Continued on next page 
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Table 4 - continued from previous page 



GAPID 


Group 


U 3 




^2x1 




Models 


3(7 


TBM 


[96, 64] 


((C 4 x C 4 ) Hp C 3 ) x v C 2 


✓ 






✓ 











[96, 65] 


Ai x v C% 


✓ 






✓ 


138 








[96, 66] 


SL(2,3) x^ C* 4 


✓ 




✓ 




48 








[96, 67] 


SL(2,3) X V C 4 


✓ 


✓ 


✓ 




48 








[96, 68] 


C 2 x ((C 4 x C 4 ) x^Gb) 


✓ 






✓ 


3648 


939 


876 


[96, 69] 


C* 4 x SL(2, 3) 


✓ 








6637 


1002 


936 


[96, 70] 


((C 2 x C 2 x C 2 x C 2 ) x v C 3 ) x,C 2 








✓ 


2433 


399 


336 


[96, 71] 


((C 4 x C 4 ) x v C 3 ) x v C 2 








✓ 


2433 


399 


336 


[96, 72] 


((C 4 x C 4 ) x v C 3 ) x v C 2 








✓ 


2433 


399 


336 


[96, 73] 


C 8 x A 4 


✓ 






✓ 











[96, 74] 


((C 8 x C 2 ) x v C 2 ) x v C 3 


✓ 


✓ 


✓ 




3884 


678 


630 


[96, 185] 


A* x v Q 8 








✓ 


16 








[96, 186] 


C 4 x 5 4 


✓ 






✓ 


113 








[96, 187] 


(C 2 x Si) x v C 2 








✓ 


16 








[96, 188] 


SL(2,3) -> C 2 ) 


✓ 








14 








[96, 189] 


C* 2 x GL(2,3) 


✓ 




✓ 




14 








[96, 190] 


(C 2 x SL(2,3)) x v C 2 










16 








[96, 191] 


SL(2,3) -^G^C 2 )^C 2 










16 








[96, 192] 


SL(2,3) -> C 2 ) x v C 2 




✓ 


✓ 




14 








[96, 193] 


(SL(2,3) x v C 2 ) x^C 2 










16 








[96, 194] 


C 2 x (A 4 x y C 4 ) 








✓ 


118 








[96, 195] 


(C 2 x C* 2 x A 4 ) xi v C 2 








✓ 


16 








[96, 196] 


C 2 x C 4 x A 4 








✓ 











[96, 197] 


D A x A 4 








✓ 


5202 


558 


486 


[96, 198] 


C* 2 x C* 2 x SL(2, 3) 










1347 


224 


198 


[96, 199] 


Q 8 x A 4 








✓ 


5187 


558 


486 


[96, 200] 


C 2 x (SL(2,3) x v C 2 ) 


✓ 




✓ 




1332 


224 


198 


[96, 201] 


(SL(2,3) ^C 2 )^C 2 










5202 


558 


486 


[96, 202] 


(C 2 x SL(2,3)) x v C 2 










5202 


558 


486 


[96, 203] 


(C 2 x C 2 x Q 8 ) x v C 3 








✓ 


189 








[96, 204] 


((Ci x D 4 ) x v ft) M„ C 8 








✓ 


189 








[96, 226] 


C 2 x C 2 x 5 4 








✓ 


42 








[96, 227] 


((C 2 x C 2 x C 2 x C 2 ) x v C 3 ) x v C 2 








✓ 











[96, 228] 


C 2 x C 2 x C 2 x A 4 








✓ 











[96,229] 


C 2 x ((C 2 x C 2 x C 2 x C 2 ) x^C 3 ) 








✓ 


4779 


853 


720 
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A. 3. Alternative Names for Some Small Groups 

In Tab. 5 we list some alternative names for the groups that we have considered in this 
publication (cf. Tab. 4). To compile this list, we have made use of refs. [10-12,21,38, 
43-45]. 

Table 5: Some aliases for the groups in Tab. 4. The first and second columns give the GAPID and the 
group names displayed by GAP, respectively. The third column shows one or more alternative names 
that are in common use in the physics and mathematics literature. The fourth column, finally, gives a 
short description of the group where appropriate. 



GAPID 


Group 


Other names 


Description 


[12,3] 




A(12),T 


Tetrahedral group 


[21,1] 


cv * v c 3 


T 7 




[24, 3] 


SL(2,3) 


V 


Double cover of A4 


[24, 12] 




A(24),0 


Octahedral group 


[24, 13] 


C 2 x Ai 


S(24) 


Pyritohedral group 


[27, 3] 


(C 3 x C 3 ) x v C 3 


A(27) 




[39, 1] 


C13 x v C 3 






[42, 2] 


C 2 x (C 7 x,C 3 ) 


Ti 4 




[48, 3] 


(C 4 x C 4 ) x v C 3 


A(48) 




[48, 28] 


SL(2, 3) -> G -> C 2 




Double cover of 5 4 


[54,8] 


((C 3 xC 3 )x v C 3 )x v C 2 


A(54) 




[57,1] 


Cl9 X ip C 3 






[60, 5] 


A 5 


S(60),7 


Icosahedral group 


[63, 3] 


C 3 x (CV x v C 3 ) 






[75,2] 


(C 5 x C 5 ) x v C 3 


A(75) 




[78, 2] 


C 2 x (C 13 x v C 3 ) 


r 2 6 




[81,7] 


(C 3 x C 3 x C 3 ) x v C 3 


£(81) 




[84, 2] 


C 4 x (CY x v C 3 ) 


r 2 8 




[93, 1] 


C 3 i x v C 3 






[96, 64] 


((C 4 x C 4 ) x v C 3 ) x v C 2 


A(96) 
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B. Construction of the Groups of Order at Most 100 



We will first describe how to generate all groups of order < 100 in GAP. Then we will 
determine which groups have a 3-dimensional irrep and/or are a subgroup of U(3) or 
SU(3). We include this information, because there seems to be a clear preference in 
model building for continuous or discrete subgroups of U(3) or SU(3). 

B.l. Generating the Groups 

The following lines of code generate the list of all groups of order < 100 using the 
SmallGroups Library [22] in GAP: 

1 SizeScreen( [ 500, ] ); 

2 groups := AllSmallGroups ( [1 . . 100] ) ; ; 

3 for g in groups do 

4 Display (StructureDescription(g) ) ; 

5 Display (IdGroup(g) ) ; 

6 chartab := Irr(g) ; ; 

7 for i in [1 .. Size (chartab) ] do 

8 Print (chartab [i] [1] , " "); 

9 od; 

10 Print("\n"); 

11 od; 

12 time; 

These lines can be entered directly at the GAP prompt. In the following we assume that 
the preceding lines have been saved in a file named smallgroups . gap that is then loaded 
and automatically executed (see line 1 below): 

1 gap> Read( "smallgroups . gap") ; 

2 1 

3 [ 1, 1 ] 

4 1 
C2 

6 [ 2, 1 ] 

7 1 1 

8 C3 

9 [ 3, 1 ] 
10 111 

We only display the first few lines of output (lines 2-10 above). For each group, there 
are 3 lines of output corresponding to lines 4, 5, 8 in the GAP code. For a non-trivial 
example, see lines 5-7 in the output. Line 5 displays the human- readable name of the 
group, line 6 gives its GAPID which uniquely identifies the group and which we will use as 
input for other GAP commands, and line 7 gives the first column of the character table, 
i.e. the dimensions of the irreps [28]. 

We find 1,048 groups of order < 100 which we list in a separate file that we have made 
available for download [23] . The first two columns of this list summarize the information 
we have obtained in this section. 
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B.2. Groups that are Subgroups of SU(3) or U(3) 



Next we determine which of these groups are subgroups of U(3) or SU(3). If a group 
is (isomorphic to) a subgroup of U(3), there is a one-to-one correspondence between 
its elements and matrices of U(3). These matrices furnish a 3-dimensional faithful rep- 
resentation of q that is not necessarily irreducible. Conversely, if q has a 3-dimensional, 
faithful representation, then q is a subgroup of U(3): For finite groups, every represen- 
tation is equivalent to a unitary representation [28], so the representation matrices are 
elements of U(3). By faithfulness, the representation p is a one-to-one mapping between 
Q and the image of p in U(3). By virtue of p being a group homomorphism, Im p inherits 
the group properties from q, and consequently Im p C U(3) is a group that is isomorphic 
to g. Finally, whether q lies in SU(3) can be verified by checking the determinant of 
representation matrices, since equivalent representations have the same determinant. 

The kernel of the representation is given by Ker p = {g E g\ char(g) = char(l)} [46], 
and thus a 3-dimensional representation p is faithful, iff 1 is the only element whose 
character is 3. For each of the 1,048 groups generated in Section B.l we calculate the 
character table. Below is the output for A 4 : 

1 gap> g := SmallGroup(12,3) ; ; 

2 gap> Display (StructureDescription(g) ) ; 

3 A4 

4 gap> chartab := Irr(g);; 

5 gap> Display (chartab) ; 

6 [ [ 1, 1, 1, 1 ] , 

7 [ 1, E(3)~2, 1, E(3) ] , 

8 [ 1, E(3) , 1, E(3)"2 ] , 

9 [ 3, 0, -1, ] ] 

On line 1, we specify the group by entering its GAPID [12,3]. Lines 6-9 give its character 
table, where E(3) = exp 27ri/3 denotes the primitive third root of unity. The first column 
gives the dimensions of the representations: 1, 1', 1", 3. On line 9 corresponding to 3, 
there is only one character equal to 3, so 3 is faithful. This proves that A4 is a subgroup 
of U(3). The representation matrices can be found by using the Repsn package [47] in 
GAP: 

1 gap> LoadPackageC'repsn") ; ; 

2 gap> for i in [1 . .Size(chartab)] do 

3 > Display (IrreducibleAff ordingRepresentation(chartab [i] )) ; 

4 > od; 

5 Pcgs([ fl, f2, f3 ])->[[[ 1 ] ] , [ [ 1 ] ], [ [ 1 ] ] ] 

6 Pcgs([ fl, 12, f3 ])->[[ [ E(3)"2 ] ], [ [ 1 ] ], [ [ 1 ] ] ] 

7 Pcgs([ fl, f2, f3 ])->[[ [ E(3) ] ], [ [ 1 ] ], [ [ 1 ] ] ] 

8 Pcgs([ fl, f2, f3 ])->[[ [ 0, 1, ], [ 0, 0,1], [ 1, 0, ] ], 

9 [ [ -1, 0, ] , [ 0, 1, ] , [ 0, 0, -1 ] ] , 
10 [ [ -1, 0, ], [ 0, -1, ], [ 0, 0, 1 ] ] ] 

Lines 8,9,10, respectively, correspond to the representation matrices of the generators 
f 1, f2, f3 of A4 for the 3-dimensional irrep. Their determinants are all 1, and thus A4 
is a subgroup of SU(3). 

In other cases, when there is no faithful, irreducible 3-dimensional representation, we 
have to consider the reducible ones. If A, B are two representations, then ohai(A(BB) = 
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char(v4) + char (5), i.e. we obtain the character of A © B by adding the rows in the 
character table that correspond to A and B. For a given group, we consider all direct 
sums that are 3- dimensional and calculate their characters. For each direct sum, the first 
element of the character will be 3, corresponding to the dimension of the representation. 
If there is more than one 3, we conclude that the direct sum is not faithful. If we cannot 
find any direct sum that is faithful, we conclude that g is not isomorphic to a subgroup 
ofU(3). 

Assume that we can find a 3-dimensional faithful, reducible representation, thereby 
proving that g is a subgroup of U(3). The representation matrices are block-diagonal, 
and each submatrix is unitary. There are two cases: All submatrices are 1 x 1, or one is 
2x2 and the other is 1 x 1. In the former case, the representation matrices are diagonal 
and commute, thus g ~ Z p x Z g x Z r C U(l) n for some n < 3. In the latter case, 
we consider the canonical embedding of the submatrices into U(3) (i.e. by extending 
the submatrix by the identity matrix to match the dimensions). Every representation 
matrix can be uniquely written as a product of these embedded submatrices, and the 
submatrices corresponding to different blocks trivially commute. This establishes that 
g is isomorphic to a subgroup of U(2) x U(l). 

We have implemented the above algorithm in a GAP script. The results have been sum- 
marized in ref. [23] and made available for download. 

B.3. Comparing Our Results to the Existing Literature 

We have compared our results to the existing literature on SU(3) subgroups [12,21,38, 
43-45]. Identifying the groups is not always straightforward, since they may appear 
under different names in different contexts, e.g. A4 is listed as A(12) in ref. [44] and as 
part of the C series in ref. [43]. 

In Tab. 1 of ref. [23] we list all groups of order at most 100 and for each group indicate 
whether it is a subgroup of U(3), SU(3), U(2), SU(2), U(2) x U(l), SU(2) x U(l), 
respectively. 

We find that the groups in our list that are subgroups of SU(3) but not of U(2) x U(l) 
agree with those in ref. [38,44,45] except in the following cases: According to ref. [38], 
the groups 05(42,2), 05(78,2), 05(84,2) are in SU(3), but our analysis along the lines of 
Section B.2 shows that they are only in U(3), in SU(3). 

Ref. [12] only explicitly lists groups that are not direct products with cyclic factors, and 
thus does not consider the groups in question, but according to Theorem II. 2 in the 
same publication, these groups are in U(3), so we have agreement. 
Also, the groups 05(36, 11), 05(72,42), 05(81,9) and 05(84, 11) were not listed in refs. [38, 
44,45], but we have verified that these groups are indeed in SU(3). 05(81, 9) and 05(84, 11) 
are part of the C series as given in ref. [43]. This has already been pointed out by ref. [12]. 
The Groups 05(36, 11) and 05(72,42) have not been explicitly listed in ref. [12], but the 
discussion following Theorem II. 2 in the same publication makes it clear that these 
groups are in SU(3). 
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Ref. [12] in Section II. 1 makes the observation that a finite subgroup of U(3) is not in 
U(2) or U(l), if and only if it has a faithful 3-dimensional irreducible representation. In 
our analysis, we find counterexamples: E.g. the group 0(16,3) C U(3) is not in U(2) 
and has no 3-dimensional irreducible representation. The reason is that the existence 
of a reducible and faithful 3-dimensional representation is already sufficient for being a 
subgroup of U(3). For more details, see Section B.2. 

For generating the A(3n 2 ) series, we have used the generators from Tab. 1 in ref. [44] 
with j = 1 and k = (also see ref. [48]). Note that if we take some arbitrary integers j 
and k, the representation may not be faithful and thus will not generate (a subgroup of 
U(3) that is isomorphic to) A(3n 2 ). Also note that in ref. [44] the generators for £(360) 
generate a group of order 1,080 which is a non-split extension 2 of Aq by C3, and not Aq. 
We agree with ref. [43] that lists the same group as £(3600). 

B.4. Groups that Contain A 4 as a Subgroup 

Since many publications in the past have highlighted A4 and its connection to tribimax- 
imal mixing, we find it useful to list the groups that contain v4 4 as a subgroup: 

1 LoadPackage (" sonata" ) ; 

2 for n in [1 . . 100] do 

3 for g in AllSmallGroups (n) do 

4 sg := Subgroups (g) ; 

5 if "A4" in List (sg,x->StructureDescription(x) ) then 

6 Print (IdGroup(g) ,"\n") ; 

7 f i; 

8 od; 

9 UnloadSmallGroupsDataO ; 
10 od; 

For every n from 1 to 100 (line 2), we generate all groups of order n (line 3). For each 
such group, we determine its subgroups (line 4) and check whether A^ is one of them 
(line 5). If the answer is positive, we print the GAPID of the respective group (line 6). 
On a technical note, since the number of subgroups becomes large, we need to increase 
the default memory allocation for GAP. The results are presented in the last column of 
Tab. 1 in ref. [23]. 



2 We are indebted to Patrick Otto Ludl for pointing out that it is not the direct product of Aq and C3 
as we had incorrectly identified in the first version of this publication due to a misinterpretation of 
the GAP output. 
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C. Breaking the Family Symmetry to Subgroups 



In general, the vev of a flavon field in an n-dimensional representation can take any value 
in IR n , but since a finite group has only finitely many subgroups, there will be a finite 
number of inequivalent vevs that break to different subgroups. It should be noted that 
the neutrino mixing angles will in general depend on the length of the vevs. Rescaling 
all vevs with the same factor, though, does not have any effect on the mixing angles. 

For definiteness, we again choose A4 x C3 as a working example, but it should be clear 
that the following discussion is completely general. Consider the GAP script listed below: 

1 SizeScreen( [ 200, ] ); 

2 LoadPackageO'repsn") ; 

3 LoadPackage (" sonata" ) ; 

4 g := SmallGroup(36, 11) ; 

5 chi:=Irr(g) ; 

6 sg:=Subgroups(g) ; 

7 Display (Size(sg) ) ; 

8 Display (List (sg, StructureDescription) ) ; 

9 for n in [1 . .Size(chi)] do 

10 PrintC'Representation" , " ",n,"\n\n"); 

11 rep : =IrreducibleAf f ordingRepresentation(chi [n] ) ; 

12 for a in sg do 

13 al : =List (a,x->TransposedMat (x~rep) ) ; 

14 Print (StructureDescription(a) , "\t\t" ,BaseFixedSpace(al) , "\n") ; 

15 od; 

16 Print (" \n\n\n" ) ; 

17 od; 

Line 4 sets the group to q = At x C3, and lines 5 and 6 calculate the character table 
and the subgroups, respectively. Line 7 prints the number of the subgroups (30 in this 
case), and line 8 lists each subgroup by its human-readable name. 
Lines 9-17 loop over the 12 irreducible representations pi of q. For each representation, 
line 11 generates the representation matrices, and lines 12-15 loop over the 30 subgroups 
t)j (j = 1,...,30) to find the vev that breaks $j to t)j as follows: For each t)j, line 
13 generates the list of the matrices M k (k = 1, . . . , dim that correspond to the 
representation pi and elements in Line 14 calculates a basis of the common eigenspace 
of the Mk for the eigenvalue 1, i.e. it calculates the vevs that are left invariant by the 
action of In the following, we will call this eigenspace V^. 

From the construction it is clear that the obtained vevs leave \)j intact, but are not 
guaranteed to break exactly to f)j, i.e. there may be another subgroup under which the 
vevs do not transform. In that case, one of the subgroups is contained in the other, 
or both are contained in a third subgroup of q. Thus, if more than one \)j leads to 
the same vev(s), the largest one is the unbroken symmetry group. This establishes a 
correspondence between the subgroups of q and the vector spaces of vevs that break 
to them, i.e. t)j -H- V\ )j (where t)j is the largest subgroup to which breaks). If there 
are no vevs that break to t)j, then VJ,. is the empty set. Note that subgroups that are 
conjugate to each other correspond to different embeddings in q and lead to different 
physics, so for our purposes they are not equivalent. To contribute to the clarity of the 
current discussion, we present the results for the irrep 3 in Tab. 6 on the following page. 
Note that V/ 10 in Tab. 6 corresponds to the vev of the flavon field </?t in Tab. 3 on page 8 
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(The flavon fields (p s and £ transform as 3' and 1", respectively, and the corresponding 
vev spaces can easily be determined by the same GAP script.) 



1)5 


= c 3 




nr 


= span { 


1,0,0), (0,1,0), (0,0,1)} 


1)19 


= c 6 






= span { 


;o,i,o)} 


1)20 


= c 6 




vf 10 

1)20 


= span { 


;o,o,i)} 


1)21 


= c 6 




K pl ° 

f)21 


= span { 


;i,o,o)} 


1)22 


= c 3 


xC 3 


K pl ° 

U22 


= span { 




1)23 


= ^3 


xC 3 


Vf 10 

U23 


= span { 


i)} 


1)24 


= c 3 


x C 3 


K pl ° 

f)24 


= span { 


:-i,i,i)> 


1)25 


= c 3 


x C 3 


vf 10 

f)25 


= span { 




1)30 


= A A 




K pl ° 

f)30 


= span { 


;o,o,o)} 



Table 6: The common eigenspaces for the eigenvalue 1 of the representation matrices for the irrep 3 that 
leave the respective subgroup i)j invariant. The irrep 3 corresponds to the 10th line of the character 
table, hence its name pio- The numbering of the subgroups t)j corresponds to their order given by the 
GAP script on p. 32. If an f)j is not listed here, this means that there exist no vevs for the irrep 3 that 
can break to this particular subgroup. 

Inspecting Tab. 6, we see another subtlety that we have to take into account (for clarity, 
we drop the superscripts). Vf, 5 is spanned by three vevs, and varying them independently 
breaks q — y f) 5 . On the other hand, Vj, 19 is a subset of Vf, 5 and breaks to a larger symmetry 
g — > hi 9 . The same is true for Vf, 20 , . . . , Vf, 25 C Vf, 5 that break to h 2 o, • • • , 1)25 ^ f)5> 
respectively. 

Assume that we choose a vev v that lies in Vf, 5 , but in none of the Vf, 19 , . . . , V(, 25 . Then 
q will necessarily break to f) 5 , since neither of the larger symmetries leave v invariant. 
Thus, we can effectively break q — y h 5 with a single vev, although Vf, 5 is 3-dimensional. 
In general, a vev will break to the symmetry that corresponds to the smallest vector 
space Vf,. in which it is contained. For two vector spaces A, B we define A < B, if 
A C B (this defines a partial ordering). A partially ordered set does not necessarily 
have a smallest element, but in this case it does. Given a vev v, there is always a 
smallest that contains it: Suppose v G fl and Vj^ (jL and <£_ . 
Then f) io ^ t)^ and v will break to a group h i2 that contains both () io and f)^ as proper 
subgroups, and as a consequence C and VJ,^ C . So either there is a smaller 
vector space in which f lies or the assumptions are not correct, i.e. Vj^ fl Vj, Si = 
or V(, io C Vt )h or Vf,^ C ^ . Thus, v is in one and only one smallest vector space. 
On a practical note, to break $j — y h, one must pick a vev v e Vf, such that f G" Kt for 
all Vm C Vf,. In the following, we will always make this assumption without explicitly 
stating it. 

In the present publication, we consider models with up to 3 flavon fields and they do 
not necessarily transform in the same representation, so in the following we extend the 
present analysis to this more general case. 
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Algorithm 1: How to find the vevs that break to a specific subgroup. 
Choose irreps p, 

,30}; 
,30}; 
30}; 



,„, p h , p i2 with io,h,i 2 e {1, . . . , 12}; 

n . , ,nv™. y t \)k for k E {1, . 



2 For p io find vev space V^ that breaks Q 



6 
7 
8 
9 
10 
11 

12 

13 
14 



3 For p ix find vev space V^ 1 that breaks Q — > t) k for A; G {1, . 

4 For pi 2 find vev space V^ 2 that breaks — >■ h^ for k e {1, . 

5 for p from 1 to 30 do 
if = then continue with next p-iteration; 
for q from 1 to 30 do 

if V^ 1 = then continue with next g-iteration; 
for r from 1 to 30 do 

if = then continue with next r-iteration; 

t = h p n t) q n fj r ; 



1 
end 



^u{(t,n;°, 



end 



15 end 



In Algorithm 1, we describe our approach to finding all inequivalent vevs that break to 
different subgroups in the case of three flavon fields. First, we choose arbitrary, but fixed 
representations p io , p in p i2 for those fields. Next, for each representation separately, we 
apply the above procedure for finding the vevs that break to all subgroups (lines 2- 
4). When all vevs are turned on, the unbroken symmetry is the intersection of h p , t) q , 
h r corresponding to the vevs VL l ° , K n , V^ 2 , respectively. Thus, for each representation 
p io , p iiy p i2 , we loop over the 30 subgroups (lines 5-15), calculate the intersections (line 
11) and collect the breaking patterns (line 12). We consider two breaking patterns 
as different, if the respective vevs Vj*° , V^ 1 , V£ do not coincide (independent of the 
subgroups t coinciding or not). 

One may wonder whether we get all the possible breaking patterns with three flavon 
fields, and the answer is positive: A single vev will necessarily breaks to a group (for a 
finite group, it is sufficient to show closure; if some elements of Q leave the vev invariant 
and are thus part of the unbroken symmetry, then this will also be true for any product 
of those elements; hence the set is closed under the group multiplication). Since q is 
finite, we can enumerate all its subgroups and find vevs that break to them (that is 
what we did in the GAP script on p. 32), thus establishing a correspondence hj <H> 
between subgroups and vevs (for simplicity, we have dropped the superscripts denoting 
the dependence on the representation). 

Suppose we are given two vevs vi, v 2 that break g — > f). Individually, v±, v 2 will break 
to some hi, h 2 , respectively, and h = hi D h 2 . Note that there may be more than one 
pair hi, h 2 whose intersection is h, and that is why we cannot consider breaking patterns 
as equivalent that lead to the same t in line 12 of Algorithm 1. The earlier established 
correspondence between subgroups and vevs now gives us hi -H- and h 2 -H- Vf, 2 , and 
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we can conclude that Vi G and t> 2 G Vf, 2 . Thus, the iteration in Algorithm 1 on the 
previous page includes the symmetry breaking pattern q — > f) for any two vevs v\ and 
V2- This concludes the proof for two vevs, and the generalization to the case of three 
vevs is straightforward. 

In this section, we have solved the problem of finding the vevs that break to all subgroups 
of a given symmetry in full detail and generality. It is important to note, though, that 
the neutrino mixing angles will not only depend on the residual family symmetry h, but 
also on the alignment and magnitude of the vevs in Vf, (as usual we assume that v does 
not lie in a smaller subspace). Resolving this ambiguity is beyond the scope of group 
theory and needs to be addressed by model building. 
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D. Clebsch-Gordan Coefficients for Finite Groups 



Currently, it is general practice to construct the Clebsch-Gordan coefficients (CGCs) for 
the various groups that are studied in physics on a case-by-case basis using heuristic 
methods. It is clear that such an approach becomes cumbersome, if one considers more 
than one group or the number of irreducible representations is large. Also, for automat- 
ing the steps from the choice of the family symmetry to finding the invariant Lagrangian 
to calculating the mixing angles and phases, we need a systematic way of deriving the 
CGCs that does not rely on the specifics of the group under consideration. 

An algorithm due to P. M. van den Broek and J. F. Cornwell [24] solves this problem 
in full generality: Given the character table and the explicit form of the unitary repre- 
sentation matrices, it calculates the CGCs for any finite group. We have implemented 
this algorithm in Python to automatically generate the CGCs for any finite group. We 
get the character table and the representation matrices from GAP that we have interfaced 
with our Python programs to achieve a high level of automization. 

To establish our notation for the CGCs and to contribute to the clarity of the discus- 
sion in Section 3, we present an explicit example of how to contract the indices in the 
tensor product of any 2 fields transforming in irreducible representations of the family 
symmetry. Consider e.g. 2 fields and if> that transform as 3's of A 4 x C3. The product 
(pi with i,j = 1, 2, 3 transforms as 

3x3=l + l'" + l (4) + 2x3 (16) 

and contains the singlet representation 1, i.e. for some choice of the product does 
not transform. The CGCs give the change of basis between the right-hand side and 
left-hand side of Eq. (16): 

0^3,3, = l k , k = l, i,j = 1,2,3 (17) 

In the general case p® q = Yl n pg x r -> the CGCs will depend on the representations p, 
q, r and on the number of times rf that r appears in the decomposition of the tensor 
product. For the case of Eq. (17), our Python script gives: 




0\ 



1 



18) 



C(p = 3,q = 3,r = l,n r pq = l)ij- k =i 

73/ 

Since 1 is a singlet, the index k only takes one value so that we can write the CGCs in 
the form of a matrix where i,j label the rows and columns, respectively. Coming back 
to the case of our 2 fields <fi and ip, we conclude that the combination 

^01^1 + 4|^ 2 + ~i ^ 3 ( 19 ) 

is invariant under A 4 x C 3 . Note that the CGCs depend on the choice of the represen- 
tation matrices. One can check the invariance of Eq. (19) by using the explicit form of 
the representation matrices obtained from the GAP script that we have discussed on p. 5. 
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We performed several checks to ascertain that the CGCs are calculated correctly. For 
one thing, we have compared our output to the (comparatively few) results that exist in 
the literature. For another, we have checked for all 90 groups in Tab. 4 on page 24 and 
for all irreducible representations p, q that the terms C i j\ z p i q^ transform as r k where 
p® q = x r - 

We find complete agreement except for A$. The problem can be traced back to the 
fact that the representation matrices for provided by GAP are not unitary. After 
choosing unitary representation matrices, the algorithm gives the correct CGCs also for 
this remaining case. 
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E. Elements of Finite Group Theory 



In this appendix we summarize some of the most important definitions and theorems 
from finite group theory that we use in the present publication. 

E.l. Direct Products 

Given 2 groups A, B we define their direct product A x B as the set of all pairs (a, b) with 
a G A and b G B, where the group operation is defined by element-wise multiplication: 

(ai, 61) • (a 2 , b 2 ) = (ai • a 2 , 61 • 62) (20) 

Conversely, when a group G is given, we can ask whether we can write it as the direct 
product of 2 of its subgroups, say A and B. From our previous definition it is clear 
that we would like 2 conditions to be fulfilled: Firstly, every element g G G should be 
expressible as a product g = a ■ b with a G A and b G B. If A and B have no common 
elements except for the identity element in G, it easily follows from the group properties 
that this decomposition is unique, and as such, we have a one-to-one correspondence 
a-b -H- (a, b). Secondly, all elements in A should commute with all elements in B so that 
we can mimic the product in Eq. (20): 

9i-92 — (ai • 61) • (a 2 • 6 2 ) = ai • bi ■ a 2 ■ b 2 = a x ■ a 2 ■ h ■ b 2 = (a 1 • a 2 ) ■ (bi • b 2 ) (21) 

With the identification a ■ b -h- (a, 6), the previous line reads: 

9i-92 = K, 61) • (a 2 , 6 2 ) = (ai • a 2 , 61 • b 2 ) (22) 

Thus we have proven that there is a one-to-one correspondence a-b -h- (a, 6) between the 
elements of G and Ax B that is compatible with the group operation, i.e. the 2 groups 
are isomorphic. For all practical purposes, we can view these 2 groups to be identical 
and write G = A x B. 

E.2. Normal Subgroups 

A subgroup N C G is called normal, if for any n G iV it holds that gng~ l G iV for all 
g G G, i.e. the operation of conjugation with an arbitrary element of G maps iV into 
itself. One then writes N <\ G. This concept is relevant in the present context, because 
both A and B are normal subgroups of G = A x B: 

ga'g- 1 = (afya'iaby 1 = aba'b^a' 1 = abb' 1 a' a' 1 = aa'a' 1 G A (23) 

The proof for B is analogous. If iV is normal in G, the cosets of G with respect to iV 
form a group, called the quotient group and denoted by G/N . Dividing by iV means 
that we identify all elements that differ by multiplication by n G N, i.e. g\ ~ g 2 , iff 
9i — n 92- In this sense, it is clear that (A x B)/B is the set of all elements of the form 
(a, 1), which is isomorphic to A, because (a, 1) ~ a • 1 = a. An analogous statement 
holds for B. 



38 



E.3. Semidirect Products 



The semidirect product is a straightforward generalization of the direct product for the 
case that A and B do not commute. Let us repeat the calculation in Eq. (21): 

9i-92 = (ai-bi)-(a 2 -b 2 ) = a r b x -a 2 -b 2 = a r a 2 -aj 1 -b 1 -a 2 -b 2 = a 1 -a 2 -(a 2 1 -b 1 -a 2 )-b 2 (24) 

Firstly, if we want to have any chance of writing g\ ■ g 2 as a product a ■ b with a £ A 
and b G B, we have to assume that a 2 1 ■ b\ • a 2 G -B . Since this must hold for all a 2 G A 
(and trivially holds for all b G B), this is equivalent to requiring that B be a normal 
subgroup. Secondly we need to know the action of A on B by conjugation which we 
denote by ip a : 

cp a :B^B, ip a (b) = a - 1 -b-a (25) 

Note that <p a : b h- )■ a -1 - 6- a is an automorphism of -B, and </? : a i— )■ </? a is a homomorphism 
from A to Aut(-B), the automorphism group of B. This may look like notational overkill, 
but it will become clear in a moment why we chose to do so. Now we can write 

91-92 = («iA) • (02,62) = («i • a 2 , ¥a 2 {bi) ■ b 2 ), (26) 

in analogy to Eq. (20), and the only difference is that the rule for multiplication gets 
slightly modified. If A and B are given as subgroups of G, the multiplication rule 
between elements from A and B is known, and we can calculate the right-hand side of 
the second term in Eq. (25). If, however, we are given 2 groups A and B that bear no 
relation to each other, we have to choose (p a from the set of automorphisms of B, and 
this serves as the definition of conjugation. In this sense, the semidirect product is not 
unique, since it depends on the choice of ip. 

To summarize, the semidirect product of A and B with respect to ip is the set of all 
pairs (a, b) with a & A and b G B, where the group operation is defined by Eq. (26) and 
ip : a i-> ip a is a homomorphism from A to Aut(-B). We use the notation by G = B x v A 
for the semidirect product. Then, (i) B < G, i.e. B is a normal subset of G, (ii) A acts 
on B by conjugation, and (iii) the quotient group G/B is isomorphic to A. Note that 
the order of the factors is significant. The semidirect product is not unique, but depends 
on the choice of <p. If f a is the identity map for all a, the semidirect product is reduced 
to the direct product. 



E.4. Short Exact Sequences 

The most general way to describe a group embedding is a short exact sequence, as we 
will now explain. An exact sequence is a collection of groups and homomorphisms 

Gl ^ G 2 ^ G3 . . . f24 Gn (27) 

such that the image of each homomorphism is equal to the kernel of the following one, 
i.e. lm(ipk) = Ker(y9fc + i). A short exact sequence is an exact sequence of the form 

1 Gl G2 Gs ^ i ; (28) 
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where 1 denotes the trivial group. A group homomorphism always maps the identity 
element onto the identity element, so Im(</?o) = 1. Because the sequence is exact, we 
have Ker(<^i) = Im(</?o) = 1, i-e. ipi is injective. Since p> 3 maps everything to 1, its 
kernel is G 3 , and by the same argument we can conclude that lm(p 2 ) = Ker(<^ 3 ) = G 3 , 
i.e. p 2 is surjective. 

The isomorphism theorem states that for any homomorphism ip : A — > B, 

lm(ip) = A/Ker(p). (29) 

It is easy to see why this holds: <p is into B, but onto Im(</?) C B, so (p : A — > lm(p) is 
surjective. The kernel Ker(tp) C A is in general not trivial, so p is not injective. Dividing 
A by the kernel (which is always a normal subset) identifies all elements in Ker(<p) with 
1, so <p : A/Ker(p) — > lm(ip) becomes injective. This concludes the heuristic proof of 
Eq. (29). 

Applied to our case, we obtain 

Im(v? 2 ) = G 2 /Ker(v9 2 ) <-> G 3 = G 2 /Im( ¥ > 1 ) ^ G 3 = G 2 jG x . (30) 

For the last equivalence we have used the fact that ipi is injective and thus establishes 
an isomorphism between G\ and Im(</?i). Again, since the kernel of a homomorphism is 
always normal, G\ ~ Im(</?i) = Ker(</? 2 ) is a normal subset of G 2 . 

We rewrite Eq. (30) using more suggestive notation: 

l^iV4G^4l A^<G and Q = G/A^ (31) 

Eq. (31) gives a description of G in terms of a normal subgroup A and the quotient 
group G/N and we say that G is an extension of Q by A. 
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